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We have explored, computationally and experimentally, the magnetic properties of (Fei-^Co^B 
alloys. Calculations provide a good agreement with experiment in terms of the saturation magne¬ 
tization and the magnetocrystalline anisotropy energy with some difficulty in describing C 02 B, for 
which it is found that both full potential effects and electron correlations treated within dynamical 
mean field theory are of importance for a correct description. The material exhibits a uniaxial 
magnetic anisotropy for a range of cobalt concentrations between x — 0.1 and x = 0.5. A sim¬ 
ple model for the temperature dependence of magnetic anisotropy suggests that the complicated 
non-monotonous temperature behaviour is mainly due to variations in the band structure as the 
exchange splitting is reduced by temperature. Using density functional theory based calculations we 
have explored the effect of substitutional doping the transition metal sublattice by the whole range 
of 5 d transition metals and found that doping by Re or W elements should significantly enhance the 
magnetocrystalline anisotropy energy. Experimentally, W doping did not succeed in enhancing the 
magnetic anisotropy due to formation of other phases. On the other hand, doping by Ir and Re was 
successful and resulted in magnetic anisotropies that are in agreement with theoretical predictions. 

In particular, doping by 2.5 at.% of Re on the Fe/Co site shows a magnetocrystalline anisotropy 
energy which is increased by 50% compared to its parent (Feo. 7000 . 3)26 compound, making this 
system interesting, for example, in the context of permanent magnet replacement materials or in 
other areas where a large magnetic anisotropy is of importance. 

PACS numbers: 75.50.Ww, 75.30.Gw, 75.50.Bb, 71.15.Nc, 71.20.Be 


I. INTRODUCTION 

Permanent magnets are used in a wide range of appli¬ 
cations, for example in data storage, energy conversion, 
magnetic refrigeration, wind power generators, to name 
just a few 52 . The determining properties of the quality of 
permanent magnet materials are the operating tempera¬ 
ture range and, primarily, the energy product. In terms 
of operating temperature range, it depends on particular 
applications, but mostly we demand that the material 
remains magnetic well above room temperature. For the 
energy product, we generally require it to be as large as 
possible. Energy product is a relatively complex quan¬ 
tity; it depends on the structure of the material (size 
of crystalline grains, their shape and orientation), but 
from the microscopic point of view it depends mainly 
on the saturation magnetization M s and magnetocrys¬ 
talline anisotropy energy (MAE), which is a measure of 
how difficult it is to rotate the magnetization direction 
by external magnetic field 1-3 . 

The best permanent magnet materials known today 
are typically utilizing rare-earth elements, most notably, 
neodymium in Nd 2 Fei 4 B and samarium in SmCos mag¬ 
nets. Both exhibit high Curie temperatures and large en¬ 


ergy products. Recently, however, the availability (and 
thus price) of the rare-earth elements became rather 
volatile, calling for development of replacement mate¬ 
rials, which would use less or none of the rare-earth 
elements. Intense research efforts have started world¬ 
wide, revisiting previously known materials, such as 
Fe 2 P 4-6 , FeNi' or Fei6N 2 8 , doing computational data 
mining among the large family of Heusler alloys 9 , ex¬ 
ploring the effects of strain 10-16 and doping by interstitial 
elements 1 ' ’ 18 , multilayers such as Fe/W-Re 19 , or as a lim¬ 
iting case of multilayers, the Llo family of compounds 20 , 
or promising Mn-based systems 21-24,50,51 , among others. 

In this work we explore the magnetic properties of the 
(Fei-^Co^B alloys. In 1970 the magnetic properties 
of these compounds have been experimentally studied 
by Iga 25 , later by Coene et al. 26 and very recently by 
Kuz’min et al. 2 '. It was found that within a certain 
range of Co contents the material exhibits uniaxial mag¬ 
netocrystalline anisotropy, with a relatively large MAE, 
significantly higher than pure iron or even cobalt ele¬ 
mental metals. To the best of our knowledge, theoretical 
understanding of the observations is missing. 

Moreover, motivated by the works of Andersson 28 and 
Bhandary 19 , we have attempted to further improve the 
MAE of the (Fei-^Co^B materials by doping them 
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with 5 d transition metal elements. The mentioned works 
have shown that the hybridization of the exchange-split 
3 d bands of Fe or Co with the 5 d, bands of heavier tran¬ 
sition metal elements may lead to a significant enhance¬ 
ment of the MAE, even if the 5 d element itself is non¬ 
magnetic or very weakly magnetic. We have thus con¬ 
sidered low percentage substitutions of Fe or Co by the 
whole series of 5 d transition metal elements and compu¬ 
tationally explored the effect of the substitutions on the 
MAE. 

Together with the first principles electronic structure 
calculations we have performed experimental studies of 
(Fei-ajCo^B alloys for selected Co concentrations and 
at the optimal concentration in terms of maximal MAE, 
we have synthesized single crystals with 2.5 at.% of se¬ 
lected 5 d elements alloyed on the 3 d site. For these single 
crystals we have then measured the magnetic properties 
and compared them to previous results 25 and our calcu¬ 
lations. 

The structure of the manuscript is the following: In 
Sec. II we present our experimental results of synthesis 
and magnetic characterization of (Fei-^Co^B alloys. 
In Sec. Ill we present our results of computational stud¬ 
ies of (Fei_ x Co x ) 2 B alloys using two different approaches 
to the problem of substitutional disorder. Different com¬ 
putational models are compared and the importance of 
correlation effects are discussied. Furthermore, we re¬ 
late the MAE to the electronic band structure around 
the Fermi energy (Ep) since it is well known that the 
MAE in d-electron systems is determined by the elec¬ 
tronic structure close to Ep 43 and in particular by the 
spin-orbit coupling (SOC) between occupied and unoc¬ 
cupied states with small difference in energy. Sec. IV 
then presents theoretical and experimental results of the 
(Fei-ajCo^B alloys doped by 5 d, transition metal ele¬ 
ments. 


II. MEASUREMENT OF MAGNETIC 
PROPERTIES 

A. Experimental Methods 

Unambiguous determination of anisotropy constants 
requires single crystals. Therefore, several single crystals 
of (Fei-^Co^B were grown as part of this work. The 
first stage consisted in melting a mixture of Fe and Co 
(both 99.99% pure) with crystalline B of 99.999% purity. 
The melting was performed in alumina crucibles placed in 
an induction furnace under a gauge pressure of 2 bars of 
pure argon. When we proceeded from the stoichiometric 
composition, the powder XRD patterns of (Fei_ x Co;c) 2 B 
alloys corresponded well to the CuAl 2 -type structure but 
the ingots also contained some precipitants (3-4 wt.%) 
rich in the 3d elements. To avoid the precipitants, the 
starting mixture was taken with a 1.5 at.% excess of B; 
the so obtained ingot appeared as single-phase in powder 
x-ray diffraction patterns and scanning electron micro¬ 


scope images. 

The ingots were subsequently remelted. For better ho¬ 
mogeneity, the melts were held for 5 minutes at 1425° C 
in an argon atmosphere. Then the temperature was re¬ 
duced slowly, at a rate of 0.1°C/min, down to 1100°C. 
At that point the power was switched off and the furnace 
quickly (initially at approx. l°C/min) cooled down to 
room temperature. The remelted ingots contained large 
(approx. 2-3 mm) crystalline grains. Individual grains 
were isolated and polished to give them the shape of 
a rectangular prism with edges parallel to the principal 
crystallographic axes. The crystals studied in more de¬ 
tail had the dimensions approx. 0.5x3x3 mm 3 with one 
of the longest edges being along [001] and the other one 
along [100]. X-ray back-scattering Laue patterns were 
used for the orientation and quality control of the crys¬ 
tals. 

The XRD measurements were done at room tempera¬ 
ture using a STOE Stadi P diffractometer with Mo Kaq 
radiation. 

Magnetization isotherms were measured in static mag¬ 
netic fields up to 3 T at temperatures ranging between 10 
and 1000 K using a Physical Property Measurement Sys¬ 
tem (PPMS-14 of Quantum Design). The magnetic field 
was applied either in the magnetization direction [001] 
or along [100]. Test measurements were also performed 
for the direction [110] and the resulting curves proved 
indistinguishable from the corresponding [100] curves. 


B. Magnetization and Magnetocrystalline 
Anisotropy 

Figure 1 displays the magnetization curves measured 
along [100] and [001] for Co 2 B (Fig. la) and Fe 2 B 
(Fig. lb) single crystals. All of M(H) curves depicted 
in Fig. 1 are shown after correction on demagnetiza¬ 
tion factor. The Co 2 B compound has uniaxial mag¬ 
netic anisotropy at low temperature (magnetic moment 
is parallel to the c-axis), however at T srt = 72 K the 
spin reorientation transition occurs, see Fig. 3. As a re¬ 
sult at T SIt the Co 2 B single crystal becomes magneti¬ 
cally isotropic and above this temperature, till the Curie 
temperature, the magnetization vector lies in the (001) 
plane (easy plane anisotropy). Contrary, the Fe 2 B single 
crystal has easy plane anisotropy at temperatures below 
T sr t = 593 K and by passing through the isotropic state 
at T srt becomes magnetically uniaxial. 

Spontaneous magnetizations, M s , determined from the 
easy axis magnetization curves are plotted versus tem¬ 
perature in Fig. 2. At low temperatures the M s was 
determined as the ordinate of the crossing-point of the 
linearly extrapolated high-field portions of the easy-axis 
curves (as shown on the 10 K curve in Fig. 1). For a more 
accurate determination of M s near the Tq (above 800 K) 
Belov-Arrott plots 29,30 were used. The continuous line 
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Figure 2: Temperature dependence of the spontaneous 
magnetization of (Fei-^Co^B. 


oblique orientation, such a field is insufficient for mak¬ 
ing the magnetization vector parallel to the applied field. 
The remaining misalignment distorts the torque curves 
and has to be corrected for. Unfortunately, there is no 
mention of such a correction in Ref. 25 . We are led to con¬ 
clude that (Fei-^Co^B are more strongly anisotropic 
than thought previously. 


Figure 1: Magnetization curves of (a) C 02 B and (b) 
Fe 2 B single crystals measured along [001] (dashed lines) 
and [100] (solid lines). 


is a fit to the following expression 2 ': 


M S {T) = M s { 0) 




(1) 

where M s (0) is spontaneous magnetization at 0 K and Tq 
is the Curie temperature. Values of the fitting parameter 
s are shown in Fig. 2. 

The anisotropy energy E a was determined as the area 
between the magnetization curves along [100] and [001] 
taken at the same temperature. No demagnetization 
correction was necessary since the sample dimensions in 
both directions had been made practically equal. The so 
obtained E a is presented as a function of temperature in 
Fig. 3. Also shown in Fig. 3 are data from an early work 
of Iga 25 . One can appreciate that Iga’s values of K 1 are 
noticeably lower than ours. The brevity of Iga’s paper 25 
precludes us from making any definite statement about 
the source of his underestimation; Iga’s A'i were deduced 
from torque curves taken in a field of 1.6 T, which is just 
above the anisotropy field at low temperatures. In an 



Figure 3: Temperature dependence of the anisotropy 
energy as determined from the leading anisotropy 
constant Ki of (Fei-^Co^B. Dashed lines are the 
values from the work of Iga 25 
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III. DENSITY FUNCTIONAL THEORY 
CALCULATIONS OF (Fex-^Co^B ALLOYS 

A. FeaB and C02B 

Theoretical analysis of (Fei-^Co^B alloys starts 
from the study of the boundary compounds Fe 2 B and 
C 02 B. Both these systems crystallize in the tetragonal 
structure with the space group I4/mcm. The calcu¬ 
lated equilibrium lattice parameters are a = 5.05 A, 
c = 4.24 A for Fe 2 B (cf. experimental values: a = 5.11 A, 
c = 4.25 A) and a = 4.97 A, c = 4.24 A for Co 2 B 
(cf. experimental values: a = 5.02 A, c = 4.22 A). 
The optimized Wyckoff positions differ only negligibly 
from rational numbers (1/6, 2/3 and 1/4). Thus, the 
calculations are based on the latter ones (see Tab. I). 
Structural optimizations are performed with the full po¬ 
tential linearized augmented plane wave method (FP- 
LAPW) implemented in the WIEN2k code 31 . The rel¬ 
ativistic effects are included within the scalar relativistic 
approach with the second variational treatment of spin- 
orbit coupling. For the exchange-correlation potential 
the generalized gradient approximation (GGA) in the 
Perdew, Burke, Ernzerhof form (PBE) 32 is used. Cal¬ 
culations are performed with a plane wave cut-off pa¬ 
rameter i?A' max = 8, total energy convergence criterion 
10 -8 Ry and with 20 x 20 x 20 k-points which was tested 
to give well converged MAE values. 

Table I: Atomic coordinates for Fe 2 B and C 02 B. 


Atom 

Site 

X 

y 


Fe/Co 

8(h) 

1/6 

2/3 

0 

B 

4(a) 

0 

0 

1/4 


The magnetocrystalline anisotropy energy (MAE) is 
evaluated as the difference between total energies cal¬ 
culated for [100] and [001] quantization axes. The ini¬ 
tial results for MAE are -0.052 meV/f.u. (-0.31 MJ/m 3 ) 
for Fe 2 B and -0.168 meV/f.u. (-1.03 MJ/m 3 ) for C 02 B. 
The experimental values of the anisotropy constant K\ 
are -0.80 MJ/m 3 for Fe 2 B and +0.10 MJ/m 3 for C 02 B, 
while K 2 are insignificant 25 . Although both theoreti¬ 
cal and experimental anisotropy constants indicate the 
in-plane anisotropy and have the same order of magni¬ 
tude for Fe 2 B, for C 02 B this comparison appears worse. 
A larger discrepancy between theoretical (-1.03 MJ/m 3 ) 
and experimental (+0.1 MJ/m 3 ) anisotropy constants 
for C 02 B occurs together with another discrepancy be¬ 
tween calculated (1.09/is) and measured (0.8/is 33 ) mag¬ 
netic moments on the Co atom. This correlation of vari¬ 
ances motivates us to study the dependence of MAE 
on magnetic moment. Fully relativistic fixed spin mo¬ 
ment (FSM) calculations were carried out based on 
the Full-Potential Local-Orbital Minimum-Basis Scheme 
(FPLO-14) 34 . These calculations were performed with a 
20 x 20 x 20 k-mesh, the PBE form 32 of the exchange- 


correlation potential and convergence criterion 10 8 Ha. 



(a) MAE (A) 



Figure 4: MAE as function of total magnetic moment 
on the 3d atom for Fe 2 B and C 02 B as calculated with 
FSM FPLO in (a). MAE(T) for Fe 2 B and Co 2 B as 
obtained by the scheme described in the text in (b). 

MAE as a function of magnetic moment (see Fig. 4a) 
indicates, for Co 2 B, that between the measured (0.8 hb) 
and calculated (1.09 ^b) values of magnetic moments, the 
MAE varies substantially in the range between +1 and 
-1 MJ/m 3 . Furthermore, the slope of MAE as a func¬ 
tion of /i at the point 1.09 [Xb is twice bigger for Co 2 B 
compared to the slope at 1.98 [is for Fe 2 B. In summary, 
FSM calculations show that the MAE of Co 2 B is a very 
sensitive function of magnetic moment on Co. Overesti¬ 
mation of the Co magnetic moment by more than 30% 
in the present calculations leads to inaccuracy of the cal¬ 
culated MAE for Co 2 B. 

Note that the dependence of MAE on magnetic mo¬ 
ment qualitatively can be related with effects of chang¬ 
ing temperature and as such it can give clues about the 
behavior of MAE as a function of temperature. If we 
assume that the magnetic moment (magnetization) de¬ 
creases with T, first slowly and then rather rapidly close 
to Tc, we may relate the calculated MAE(/j) to experi¬ 
mental MAE(T). For Fe 2 B we see that with decreasing 
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of fj, the MAE starts from a negative value then increases 
to a positive maximum and after that it decreases to¬ 
wards zero, in a qualitative agreement with the exper¬ 
imental Ki(T) curve. For Co 2 B, if we start from the 
experimental [i value, we see that MAE starts from a 
positive value, decreases to negative minimum, and then 
once again increases, also in qualitative agreement with 
experimental K\(T). The analysis above should however 
be taken with some care since the curves in Fig. 4a are 
from collinear configurations in which each atomic mo¬ 
ment decreases, while in the experiments non-collinear 
configurations certainly also play role. Still, the observa¬ 
tion described above gives an indication to the source of 
the rather complicated temperature dependence of the 
MAE in this system. It appears that it might be ex¬ 
plained by the variation in the MAE due to changes in 
the band structure as the exchange splitting is reduced 
with increasing temperature. 

In an attempt to quantify the above arguments we 
present a highly simplified model to map the MAE(/x) 
curves in Fig. 4a to temperature. To large extent the 
temperature dependence of the MAE is expected to come 
from two sources, namely the reduction in the spin split¬ 
ting of the bands and the fluctuations in the directions of 
the moments. The magnetism observed here is expected 
to be of a highly itinerant character and since the major¬ 
ity spin channel contains unoccupied states, it is itinerant 
ferromagnetism of the weak kind. Such a system is ex¬ 
pected to exhibit a magnetic moment which scales with 
temperature approximately as 4, 

rri 2 

( 2 ) 

1 c 

Experimental values for Curie temperatures allow an ini¬ 
tial mapping of MAE to temperature from the data in 
Fig. 4a via Eq. 2. In the next step we wish to take 
into account also the directional fluctuations of the spins. 
The Callen and Callen model 48 predics that such fluctu¬ 
ations result in a (m(T)/m(0)) 3 behaviour in the MAE. 
A classical mean field model with a Langevin function 
is used to enforce such a scaling. Finally, we choose the 
the zero-temperature reference values for the moments 
in order to obtain MAE(T = 0) values closer to experi¬ 
ment and hence a quantitatively somewhat more accurate 
model than that obtained by using the experimentally or 
computationally obtained moments at low temperature. 
In doing so the zero-temperature moments are chosen 
to mco(T = 0) = 0.7/tb and mp e (T = 0) = 2.1/zb- 
The resulting MAE(T) curves for Fe 2 B and C 02 B are 
included in Fig. 4b. It appears that this crude and sim¬ 
plified model at least captures the main features of the 
experimental MAE(T) curves when comparing to Fig. 3. 
This further supports the idea that much of the interest¬ 
ing non-monotonous behaviour in the MAE(T) is due to 
variations in the band structure as temperature reduces 
the exchange splitting. 


B. Virtual Crystal Approximation 

To study the whole range of compositions between 
Fe 2 B and C 02 B the Virtual Crystal Approximation 
(VCA) method was used. VCA imitates a random on¬ 
site occupation of two types of atoms by using a vir¬ 
tual atom with an averaged value of the nuclear charge 
ZycA- In the case of the (Fei-^Co^B alloys it is 
Zy ca = (1 — x) ■ Zp e + x ■ Zq 0 = 26 + x. The num¬ 
ber of valence electrons on the VCA site is also increased 
correspondingly. 



Figure 5: Lattice parameters as functions of x in 
(Fei-^Co^B, from experiment and calculated with 
WIEN2k in the GGA, treating disorder by the VCA. 


In the first step of investigating the alloy compound, 
the lattice parameters are studied as a function of alloy 
concentration using WIEN2k within the GGA. Calcula¬ 
tions were performed with and without the inclusion of 
SOC which, however, had a negligible effect on lattice 
parameters. Fig. 5 shows experimentally measured and 
theoretically calculated lattice parameters, a and c, as 
functions of x in (Fei_ x Co x ) 2 B. 

Overall, the discrepancy between theory and experi¬ 
ment is small and in general the variation in lattice pa¬ 
rameters are small, in the order of 1-2% over the whole 
range. Calculations were also performed within the lo¬ 
cal density approximation (LDA) 36 with similar results 
for a and c but somewhat greater disagreement with ex¬ 
periment (not shown). Further results presented in this 
section are obtained with structures from linear interpo¬ 
lation between theoretically obtained lattice parameters 
of the end compounds. For a few alloy concentrations 
the MAE obtained by using the computationally opti¬ 
mized lattice parameters was compared to that obtained 
by using the lattice parameters from linear interpolation 
and the difference was found to be negligible. It was also 
investigated how internal atomic positions vary with x 
and the variation was found to be minute, within numer¬ 
ical and experimental uncertainties, whereby fixed values 
were assumed as presented in Tab. I. 

Fig. 6 reveals an excellent agreement between MAE’s 
as functions of x calculated with the WIEN2k and FPLO 
codes. Such consistency reflects a good convergence and 
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Figure 6: MAE as a function of x in (Fei-^Co^B 
calculated with WIEN2k and FPLO, in both cases 
treating disorder by the VC A. 


high numerical accuracy of these results. The func¬ 
tion MAE(a;) is in good qualitative agreement with mea¬ 
sured Ki(x) at low temperatures 25 . Both the uniaxial 
anisotropy maximum, closer to Fe-rich side, and the in¬ 
plane anisotropy minimum, on Co-rich side, are repro¬ 
duced. The highest uniaxial anisotropy, with MAE = 1.5 
MJ/m 3 , occurs for the (Feo.eCoo.^B alloy. Quantita¬ 
tive overestimation of the MAE for the whole concentra¬ 
tion range brings reminiscence to the case of tetragonally 
strained Fe/Co alloys, where VCA treatment of substi¬ 
tutional disorder 10 also led to an overestimation of MAE 
compared to experiments 12 or more sophisticated treat¬ 
ments of disorder 15,16 . Therefore, in the next subsection 
we will also present results based on a more advanced 
description of disorder, the coherent potential approxi¬ 
mation (CPA). 

We continue first with our analysis of the VCA re¬ 
sults by exploring the relation of the MAE to the orbital 
angular momentum anisotropy. Calculated orbital mag¬ 
netic moments hl vary between 0.008 and 0.014 /rs per 
Fe/Co atom. As seen from Figs. 6 and 7 the MAE closely 
follows the orbital moments anisotropy, as suggested by 
Bruno’s formula 35 , although with a slight shift towards 
the positive values. 

Total magnetic moment hs+l per (Fe/Co) virtual 
atom calculated as a function of x in (Fei_ x Co a; ) 2 B 
is compared with experimental values of mean 
magnetization 33 in Fig. 8. Discrepancy between experi¬ 
ment and theory is equal to 0.1 ^ib /atom for Fe 2 B and 
0.3 [IbI atom for C 02 B. For intermediate concentrations 
discrepancies scale almost linearly with x. 

In Fig. 9 the FSM and VCA dependencies of MAE 
are presented as a 2D color map in superposition with 
the theoretical equilibrium magnetic moments ^ s+l s de¬ 
noted by black dots. The MAE landscape reveals that 
by going from magnetic moment of pure Fe 2 B towards 



Co concentration x 



Co concentration x 

Figure 7: The orbital magnetic moment /i£ and 
difference of ^l’s for [100] and [001] quantization axis 
as functions of x in (Fei-^Co^B, calculated with 
FPLO treating disorder by the VCA. 


the much lower magnetic moment of C 02 B the MAE 
path passes across a broad range of positive values and 
a narrower and steeper valley of negative values. Par¬ 
ticularly, the final value of MAE on the C 02 B side de¬ 
pends sensitively on the Co magnetic moment. Clearly, 
on the Co-rich side the overestimation of Co moment af¬ 
fects substantially the shape of MAE curve. Focusing on 
the uniaxial anisotropy, the highest MAE value along the 
equilibrium path is reached for the (Feo.eCoo.^B alloy. 

We note that on the calculated MAE landscape, there 
is a region with about three times higher values of MAE 
for /jl about 1.25 ^s/atom on the Fe-rich side. To ap¬ 
proach this region, starting, e.g., from the (Feo.eCoo.^B 
composition, both the magnetic moment and the d- 
electrons number (proportional to x) should be reduced. 
The desired magnetic moment reduction from 1.75 to 
around 1.25 ^s/atom is around 25%. This can be ob¬ 
tained by alloying Feo. 6 Coo .4 with 25% of a suitable non- 
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Co concentration x 

Figure 8: Total and species-resolved magnetic moments 
as functions of x in (Fei-^Co^B. Experimental values 
from Ref. 33 , theoretical moments calculated with FPLO 
treating disorder by the VC A. 


magnetic element. The non-magnetic alloying should at 
the same time decrease x (which may be understood as 
the number of d electrons beyond those of Fe) by about 
0.2. Such decrease can be obtained by alloying Feo. 6 Coo .4 
with elements having less of d electrons, i.e., the elements 
from columns of the periodic table preceding the column 
with Fe. Among 3d elements which fulfill the latter con¬ 
dition are Cr or Mn, but these often carry rather large 
magnetic moments, therefore they would not fulfill the 
first condition - reduction of average moment per atom. 
From 4 d and 5 d elements, possible candidates are Mo 
and Tc, or W and Re, respectively. Doping with heavier 
transition metal elements will be explored in Sec. IV. 


C. Coherent Potential Approximation 

As seen in the previous section, the VCA yields a cor¬ 
rect qualitative agreement for MAE as a function of x in 
(Fei_ a; Co 2 .) 2 B when comparing with experiments. How¬ 
ever, as previously discussed, it quantitatively overes¬ 
timates the MAE. Hence, calculations have also been 
performed using the coherent potential approximation 
(CPA) to treat the alloying. These calculations were car¬ 
ried out with the SPRKKR 39,40 method and GGA 32 for 
the exchange-correlation potential. The MAE was eval¬ 
uated by the torque method 41 and as much as 160000 
fc-vectors were used for the numerical integration over 
the full Brillouin zone in order to obtain MAE values 
with numerical accuracy within a few percent. Calcula¬ 
tions were done using the lowest temperature, i.e. 77 K, 
experimental lattice parameters reported by Iga 25 . 

Fig. 10 shows the MAE(i) and saturation magneti¬ 
zation, M s (x), as result of these calculations. As ex¬ 
pected, application of CPA instead of VCA removes the 


MAE [MJ/m 3 ] 



Figure 9: Fully relativistic FPLO calculations of MAE 
as a function of x and total magnetic moment (/ jls+l) 
on 3d atom for (Fei-xCo^B. Disorder was treated by 
the VCA and hs+l were stabilized with fixed spin 
moment FSM approach. Equilibrium hs+l ’s (see 
Fig. 8) are denoted by black dots. 


problem of significantly overestimating the MAE of the 
alloy, at least on the Fe-rich side of the diagram. In¬ 
stead the qualitative behavior appears to be inconsistent 
with the experiments of Iga 25 , as well as the experimen¬ 
tal data provided in this paper, on the Co-rich side of 
the diagrams. This failure is most likely due to the 
atomic sphere approximation (ASA) which approximates 
the potential as being spherically symmetric around each 
atom, in contrast to the VCA calculations presented in 
the previous section, which included full potential (FP) 
effects. This is in contrast with a recent publication 42 , 
where authors suggest that the qualitatively wrong be¬ 
havior in the MAE predicted for the Co-rich side is a 
failure of exchange-correlation functional itself, yielding 
too high magnetization. Nevertheless, the CPA descrip¬ 
tion in ASA coincides with experiments on the Fe-rich 
side. In particular the region around 0.2 < x < 0.4, with 
rather large uniaxial MAE, which is the most interest¬ 
ing region from a practical perspective, is correctly de- 











scribed. The maximum MAE is found for x = 0.3 where 
MAE = 131 fieV /f.u. = 0.77 MJ/m 3 . 




Figure 10: MAE and saturation magnetization as 
functions of x in (Fei-^Co^B as calculated with 
SPRKKR treating disorder by the CPA. 


In order to further investigate the importance of full 
potential (FP) effects, calculations using both ASA and 
FP in SPRKKR were performed for Fe 2 B and C 02 B. The 
band structures around the Fermi energy resulting from 
these calculations are shown in Fig. 11. In the case of 
Fe 2 B, there are quantitative changes in the bands but 
qualitatively they look similar. The MAE of the system 
changes from —0.11 meV/f.u. = —0.64 MJ/m 3 in ASA 
to —0.22 meV/f.u. = —1.28 MJ/m 3 in FP, which can 
be compared to the experimental value of —0.80 MJ/m 3 . 
This variation in MAE is certainly of numerical signifi¬ 
cance but still the numbers qualitatively agree on the sign 
and order of magnitude and both are reasonably close to 
experiment. 

For the case of C 02 B, the MAE changes from 
—0.89 meV/f.u. = —5.4 MJ/m 3 to —0.37 meV/f.u. = 
—2.2 MJ/m 3 upon inclusion of FP in the calculations. 
The later value is significantly closer to the experimen¬ 


tal low temperature values, as well as to the results of 
FP calculations presented in previous sections of this pa¬ 
per and it could potentially allow for reproducing the 
correct qualitative behavior for the Co-rich alloys. In 
the band structure one can observe certain qualitative 
changes which can be of relevance with respect to the 
MAE. In particular around the high symmetry point X, a 
relatively flat band which is unoccupied in ASA becomes 
occupied in FP which is relevant to the MAE since bands 
crossing the Fermi energy tend to dramatically affect the 
MAE 44 . However, to allow for further analysis of the 
importance of various bands to the MAE, we have used 
scalar relativistic, spin polarized calculations, neglecting 
spin-orbit interaction, carried out in WIEN2k to identify 
the orbital and spin characters of some bands which are 
marked in Fig. lib. This information can be used in an¬ 
alyzing how the coupling between the bands contribute 
to the MAE by considering the SOC in second order per¬ 
turbation theory 35 , where the only non-zero contribution 
comes from coupling between occupied and unoccupied 
states, and by looking at the spin-orbit matrix elements 
which are tabulated e.g. in Ref. 45 . In the bottom part 
of Fig. lib, the WIEN2k FP-LAPW bands, which are 
very similar to the SPRKKR-FP bands, are also shown 
for comparison together with the MAE contribution per 
k-point as obtained by the magnetic force theorem 44,46 . 
In the ASA case, the highest occupied band at X shows 
mainly minority spin d z 2 character while the lowest un¬ 
occupied is minority spin d xz . These states do not couple 
via spin-orbit interaction whereby a weak contribution to 
the MAE is expected in this region. As the minority spin 
d xz band becomes occupied in FP, it will couple to both 
of the unoccupied states with magnetic quantum number 
m = ±2 which have opposite spin character and cause 
MAE contribution of opposite signs resulting in cancel¬ 
lation and again a weak effect on the MAE as also seen 
explicitly in the bottom panel of Fig. lib. An important 
negative MAE contribution instead appears to come from 
a region around T and the coupling of the occupied mi¬ 
nority spin d xz band with the flat unoccupied minority 
spin d xy band (in general, coupling between states with 
the same spin but different magnetic quantum numbers 
contribute negatively to the MAE 44,45 ). These two bands 
are closer to each other in the ASA case allowing for this 
negative region to be overestimated and contribute to the 
too large negative MAE obtained for C 02 B in ASA calcu¬ 
lations. Typically, flat bands such as that just above the 
Fermi energy at T can contribute strongly to the MAE 
and small shifts in such bands will be essential and might 
partially explain the sensitivity in the MAE of C 02 B to 
various parameters. Including FP effects, the MAE is 
still more strongly negative both than the experimental 
value, which is slightly positive at low temperatures, and 
that from previous FP (WIEN2k and FPLO) calcula¬ 
tions, which yield a negative value with approximately 
half the magnitude. However, the previous calculations 
were performed using computationally optimized lattice 
parameters while the SPRKKR calculations used slightly 
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Figure 11: Band structure of Fe 2 B and C 02 B obtained 
in SPRKKR-ASA and SPRKKR-FP and for Co 2 B also 
that obtained with FP-LAPW calculations, in the lower 
panel,together with the MAE contribution of each 
k-point (solid red line) as obtained by the force 
theorem. The orbital and spin characters of some bands 
are indicated with majority and minority spin denoted 
by t and J,, respectively. 


larger experimental lattice parameters. Calculating the 
MAE for C 02 B with the computationally optimized lat¬ 


tice parameters presented in previous sections of this pa¬ 
per in SPRKKR-FP yields MAE = -0.19 meV/f.u. = 
— 1.2 MJ/m 3 , in good agreement with the FPLO cal¬ 
culations (—0.17 meV/f.u. = —1.03 MJ/m 3 ). Another 
indication that FP effects are of great significance for 
C 02 B comes from the density of states at the Fermi en¬ 
ergy. For C 02 B this quantity changes significantly from 
8.1 to 5.9 states/eV upon inclusion of FP effects, while for 
Fe 2 B the change is from 4.4 to 4.3 states/eV. Moreover, 
the relatively large density of states at the Fermi energy 
for C 02 B could partially explain why the MAE of this 
compound appears more sensitive to various parameters 
since, with many states near the Fermi energy, a precise 
description of these states will be essential for a correct 
MAE. In conclusion it is clear that the MAE, which is in 
general a challenging quantity to obtain from first prin¬ 
ciples, is for this system, especially on the Co-rich side, 
particularly sensitive to various parameters and difficult 
to describe correctly. However, the Fe-rich side of the 
compound appears to be well described using ASA and 
CPA, which will hence be used in Sec. IV. Since the full 
potential VCA calculations overestimate the MAE but 
produce the qualitatively correct shape of the MAE ( 2 ;) 
curve, while the ASA CPA calculations yield a qualita¬ 
tively incorrect curve on the Co-rich side but the correct 
behaviour on the Fe-rich side, one is led to believe that 
a combination of FP and CPA might produce an MAE 
curve which is both qualitatively and quantitatively in 
reasonable agreement with experiment for all alloy con¬ 
centrations. 

For further insight into the MAE of the alloy system we 
have plotted the Bloch spectral functions from SPRKKR- 
CPA calculations together with the band structure from 
FPLO-VCA calculations for a few alloy concentrations 
in Fig. 12. In addition to the smearing in the spec¬ 
tral functions, which can be seen to increase somewhat 
with x , certain bands appear slightly shifted in compar¬ 
ing the bands and spectral functions. These shifts are 
expected to be mainly due to FP effects, which are ab¬ 
sent in SPRKR-CPA, a conclusion which is supported 
by comparison of Fig. 12 to Fig. 11a, and in general 
the electronic structure appears very similar in the two 
models. This similarity is the reason why the VCA and 
CPA tend to yield qualitatively similar results while the 
quantitative differences appear because the band smear¬ 
ing is absent in the VCA. In Ref. 42 the filling of states 
with alloy concentration was investigated and related to 
the variation of the MAE. It was suggested 42 that the 
spin-diagonal part of the spin-orbit operator is most im¬ 
portant with rather strong and nearly constant, positive 
MAE contributions coming from the coupling between 
majority spin states while coupling between minority spin 
states yields the variation in the MAE. More specifically, 
it was suggested that minority spin bands yield a nega¬ 
tive contribution for x < 0.3, which is suppressed as the 
two minority spin bands on opposite sides of Ep at the T 
point for x = 0.2 become occupied, which they both are 
at x = 0.3. Identification of the character of the bands in 
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(a) x = 0.1 


(b) x = 0.2 


(c) x = 0.3 


Figure 12: Spectral functions from SPRKKR in the CPA and electronic bands from FPLO in the VCA for 
(Fei-^Co^B. In (a) the main spin and d-orbital character of the two lowest unoccupied and two highest occupied 
states at T are indicated with 4- denoting the minority spin channel. 


Fe 2 B using a scalar relativistic, spin polarized calculation 
neglecting SOC supports that the two relevant lowest un¬ 
occupied bands at T in Fig. 12a are mainly of minority 
spin character with magnetic quantum number m = ±1. 
It is also confirmed that the two highest occupied states 
at r are of minority spin character with magnetic quan¬ 
tum number m = ±2. Since the spin-diagonal part of the 
spin-orbit coupling between states of different magnetic 
quantum numbers yiels a negative MAE contribution, 
this coincides with the picture given in Ref. 42 . How¬ 
ever, for a more complete picture of the total effect on 
the MAE, contributions from bands around the Fermi 
energy in the entire Brillouin zone should be carefully 
considered. This is illustrated further by the observation 
that the qualitative description, just discussed to explain 
why there is a maximum in the MAE for x = 0.3, appears 
to hold also for the FPLO VCA calculations even though 
these do not show the maximum MAE until x = 0.4. We 
finalize this discussion by concluding that the MAE is a 
delicate property sensitive to the electronic structure in 
a region near the Fermi energy while major changes tend 
to be traceable to bands crossing the Fermi energy. 


D. Electron Correlations in Fe2B and C02B within 
DMFT 


The disagreement between the presented DFT calcula¬ 
tions and experimental measurements regarding the mag¬ 
netic moment of C 02 B motivates further investigation to 
resolve this issue and correlation effects beyond standard 
DFT methods will therefore be considered. In order to in¬ 
vestigate correlation effects on the magnetic properties of 
the Fe 2 B and C 02 B compounds, we have performed elec¬ 
tronic structure calculations within dynamical mean-field 
theory (DMFT), 55,56 using the so-called DFT+DMFT 
approach 53,54,57,58 as implemented in the full-potential 
linear muffin-tin orbital (FP-LMTO) code “RSPt”. 59,60 
The 3d subset of correlated orbitals was defined using 
the “muffin-tin head” (MT) projection. For further de- 



Figure 13: Projected density of states of the 3d states 
in Fe 2 B (fig. a) and C 02 B (fig. b), respectively, 
calculated within DFT and DMFT (SPTF solver). The 
Fermi level is set to zero. 


tails we refer the reader to Refs. 58, 61, and 62. The 
DMFT calculations are performed at a finite temperature 
set to room temperature in this study. Spin-orbit cou¬ 
pling in the muffin-tin region was included in the calcu¬ 
lations. The impurity problem occurring within DMFT 
was solved via the spin-polarized T-matrix fluctuation- 
exchange (SPTF) solver, 63,64 which is based on pertur¬ 
bation theory. The Hubbard U was set to 1.5 eV and 
2.5 eV for Fe 2 B and Co 2 B, respectively, and the Hund’s 
exchange J to 0.9 eV. These values of U are chosen since 
values in the range 1 — 3 eV are commonly used for metal¬ 
lic Fe and Co while using a somewhat larger value for Co 
than for Fe has been shown to yield a good description of 
the orbital magnetism in these elements 49 . The around 
mean-field approximation was used for the double count¬ 
ing correction. 

The calculated magnetic moments within DFT and 
DMFT are given in Table II. DMFT brings a decrease in 
the spin polarization for the transition metal moments in 
both Co 2 B and Fe 2 B. The decrease in the Co spin polar¬ 
ization is considerable, in good agreement with the exper- 
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imentally measured moments. The decrease can be un¬ 
derstood by investigating the projected density of states 
of the transition metals, shown in Fig. 13. Negligible 
differences among DFT and DMFT are found for the mi¬ 
nority spins for both compounds, as well as the majority 
spins for Fe 2 B. As for the C 02 B majority spins, the peak 
below the Fermi level is pushed to higher energies, thus 
decreasing the overall spin polarization. In addition, a 
satellite-like feature is observed at low energies. 

Based on these results we conclude that not only an ac¬ 
curate description of the electrostatic potential, but also 
a treatment of correlation effects more advanced than 
that offered by the GGA is of importance to obtain a 
satisfying picture of the magnetic properties of C 02 B. 
For Fe 2 B on the other hand, such effects appear to be of 
less significance. 


IV. SUBSTITUTIONAL DOPING BY 5 d 
TRANSITION ELEMENTS 

Various ways of enhancing magnetocrystalline 
anisotropy in 3 d magnets have been discussed by 
Kuz’min et al. 2 ‘. In this work we tried to prove the 
concept of 3d-5d interactions in a more complex system 
like (Fe 0 . 7000 . 3 ) 26 . As spin-orbit coupling is essential 
for the MAE, heavier elements with large SOC can 
possibly have significant effects on the MAE. Even 
non-magnetic elements can be of great importance via 
hybridization 19,28 . Hence, the possibility of tailoring 
the MAE by adding small amounts of 5 d elements 
substituting some of the Fe and Co in the concentration 
around (Feo. 7000 . 3 ) 26 , where the MAE is large and 
positive, has been explored. 


A. Theory 

Using the same computational methods as described 
in Sec. IIIC (SPRKKR-ASA in the CPA), the magnetic 
properties and in particular the MAE, was calculated 
for various elements X from the 5 d row of the peri¬ 
odic table in the compounds (Feo. 675000 . 3 X 0 . 025)26 and 
(Feo. 675000 . 275 ^ 0 . 05 ) 26 . The result is shown in Fig. 14 
with the X marked on the .T-axis and a dotted line indi- 


Table II: Calculated magnetic moments within DFT 
and DMFT (SPTF solver). 


Co 2 B 

Co spin moment [/ib] 

Co orbital moment [/xb] 

DFT 

1.083 

0.036 

DMFT 

0.813 

0.034 

Fe 2 B 

Fe spin moment [/xb] 

Fe orbital moment [/xb] 

DFT 

1.938 

0.041 

DMFT 

1.841 

0.041 


eating the MAE of (Feo. 7000 . 3)26 for comparison. The 
calculations indicate that doping with W or Re indeed 
appears to cause a remarkable increase of the MAE, 
with the cost of a small decrease in M s due to substi¬ 
tuting some magnetic elements for non-magnetic ones. 
The greatest increase in MAE is for Re where there is a 
rise from MAE = 0.77 MJ/m 3 to MAE = 1.14 MJ/m 3 
or 1.38 MJ/m 3 , assuming 2.5 % or 5.0 % substitutions 
per Fe/Co atom, respectively. 



Figure 14: MAE for various elements X in 

(Feo. 675000 . 3 X 0 . 025)26 and (Feo. 675000 . 275 ^ 0 . 05 ) 26 . 
The dotted line indicates the MAE of (Feo. 7000 . 3)26 for 
comparison. Crosses indicate the MAE obtained when 
the spin-orbit coupling on the 5 d elements is set to zero. 

The large increase in MAE only appears to be seen 
when doping with W or Re and not with other elements 
studied, in agreement with qualitative prediction based 
on Fig. 9. The MAE is sensitive to the electronic struc¬ 
ture around the Fermi energy and the dopants will not 
only contribute with stronger SOC but also affect the 
general electronic structure. In some cases this effect will 
be beneficial for the MAE, such as for W or Re here, and 
in some cases not, such as for Ir or Ag where the MAE 
is reduced. In order to investigate the importance of the 
strong SOC of the 5 d atoms, calculations were performed 
where the SOC is set to zero on these atoms and the re¬ 
sults are included in Fig. 14. There is a small variation in 
the MAE also in this case where the SOC is zero on the 
5 d atoms and the trend in this variation is essentially the 
same as in the case with all SOC included. However, it is 
clear that these variations can be significantly enhanced 
due to the strong 5 d SOC and that this is what yields 
the immense MAE for Re or W doping. Other relevant 
changes to the electronic structure around the Fermi en¬ 
ergy caused by the dopants could in principle be analyzed 
in the spectral functions but it is difficult to observe spe¬ 
cific changes due to the small amounts of dopants in the 
rather complicated spectral functions for the disordered 
systems and hence these have not been included. As the 
effect of the dopants is more complex than simply caus- 




















12 


ing an enhancement of the MAE, it could be of interest 
to study 5 d substitution around other values of x than 
x = 0.3 as well. 


B. Experiments on (Fe 0 .675Co 0 .3Xo.o25)2B 

As predicted by the theoretical calculations we tried 
to substitute the pure (Feo.rCoo.a^B system with 5 d el¬ 
ements of W, Re and Ir. The room temperature XRD 
measurements of 5 d doped samples are shown in Fig. 15. 



29 (deg.) 

Figure 15: Room temperature XRD plots of 
(Feo. 675 Coo. 3 Xo. 025)22 for X = (a) W, (b) Ir and (c) 
Re. Colored lines indicate the Bragg positions of 
Fe 2 B-phase and CoWB-phase in red and blue, 
respectively. The inset in (a) shows an enlarged region 
of 17° — 23° for W substituted sample. 

The room temperature XRD pattern (Fig. 15a) of 
W doped sample shows two different phases, main 


(Fe,Co) 2 B phase and a secondary CoWB (Pnma) phase. 
We tried different synthesis methods to avoid the forma¬ 
tion of the secondary phase, but in all cases W tends to 
create a secondary phase rather than replacing the Fe in 
the main phase. In contrary to the W substitution, single 
phase materials were obtained for Ir (Fig. 15b) and Re 
(Fig. 15c) substitutions. 



Figure 16: Magnetization curves of (Feo. 7 Coo. 3 ) 2 B and 
(Feo. 67 . 5 Coo. 3 X 0 . 025)22 single crystals along [100] and 
[001] (X = Re, Ir). 


The comparison of the room temperature magnetic 
measurements of pure (Feo. 7000 . 3)22 and 5 d substituted 
(Feo. 675000 . 3 X 0 . 025)22 with X = Re, Ir samples are 
shown in Fig. 16. For the Re substituted sample, the 
saturation magnetization at 2 T decreases from 143.3 
Am 2 kg _1 to 122.2 Anrkg -1 which is expected due to 
reduced amount of magnetic elements and, consequently, 
increased Fe-Fe interaction distances. Contrary to the 
observed decrease in the saturation magnetization, quite 
a big improvement is observed in the anisotropy field. 
The Re substitution increased the anisotropy field of the 
pure system from 1 T to 1.6 T. 

The Ir substitution slightly decreases the saturation 
magnetization at 2 T to 138.2 Am 2 kg _1 (Fig. 16). In 
addition the anisotropy field is also decreased down to 
0.8 T, which leads to a decrease in the leading anisotropy 
constant. 

MAE values for both pure and 5 d substituted samples 
are shown in Fig. 17 for different temperatures. The re¬ 
sults show that for the measured temperature interval an 
increase in MAE is observed for the Re substitution. The 
room temperature MAE increases from 0.42 MJm -3 to 
0.63 MJm -3 for Re substitution (approx. 32% increase). 
This ratio stays almost constant down to 10 K. This rel¬ 
ative increase is similar to the theoretical data shown in 
Fig. 14. 

The MAE for the Ir substituted sample is 0.36 MJm ~ 3 
at room temperature, which is approximately 15% lower 
than for the pure system. For the high temperature re- 
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Figure 17: Temperature dependence of the leading 
anisotropy constants of (Feo. 7000 . 3 ) 26 , 
(Feo. 675 Coo. 3 R- e o. 025)26 and (Feo. 675000 . 360 . 025 ) 26 . 


gion (T > 600 K), On the other hand, the MAE shows 
higher values compared to the pure system. 


V. CONCLUSIONS 

In summary, we have presented a combined theoreti¬ 
cal and experimental study of structural and magnetic 
properties of (Fei-zCo^B alloys with emphasis on the 
magnetocrystalline anisotropy as this property is essen¬ 
tial, for example, for replacement material candidates 
for permanent magnet applications. The qualitative 
shape of the MAE curve as function of Co-concentration 
agrees with experiment for full-potential calculations in 
the VCA. However, the quantitative discrepancy is se¬ 
vere as commonly observed when studying the MAE in 
the VCA 15-1 '. Calculations in the atomic sphere approx¬ 
imation, treating disorder in the CPA, solves the prob¬ 
lem of the quantitative disagreement with experiment on 


the Fe-rich side of the alloy but instead yields a quali¬ 
tatively incorrect behavior on the Co-rich side. Further¬ 
more, DFT with the generalized gradient approximation 
for exchange-correlation potential results in an overesti¬ 
mate, compared to experiment, of the magnetic moment 
of C 02 B, which is shown to strongly affect the MAE. 
It is shown that a reduction of the exchange splitting 
and the correct magnetic moment of C 02 B is reproduced 
if electron correlations are treated by dynamical mean 
field theory. This leads us to believe that it would in 
principle be possible to correctly describe both the mag¬ 
netic moment and the MAE over the whole range of alloy 
concentrations by taking into account both full potential 
effects and dynamical mean field theory while describing 
the alloy with the coherent potential approximation. Un¬ 
fortunately, such calculations are outside of our current 
capabilities. This provides an alternative description of 
the observed MAE curve compared to a recent work 12 , 
where authors artificially reduce the exchange interaction 
in order to describe the upturn in MAE on the Co-rich 
side. In addition, we have studied the effect of doping 
these materials with 5 d elements on their magnetic prop¬ 
erties. We have found this to be a viable route to en¬ 
hancing the MAE of (Fei-^Co^B alloys. Simulations 
suggest that 5 at.% doping by W or Re should approx¬ 
imately double their MAE, with only modest reduction 
of saturation magnetization. Experiments support this 
finding for Re, while synthesis of W-doped (Fei-a,Co^B 
was not successful due to formation of additional phases. 
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